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Abstract 
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1 Introduction 

A fundamental feature of the Standard Model 
and its three quark generations is that all CP 
violation phenomena are the result of a single 
phase in the CKM quark- mixing matrix pQ. The 
validity of this model may be tested in several 
ways, and one — verifying the unitarity condi- 
tion V ud V* b + V cd V* b + V td V* b = — is readily 
applicable to B mesons. This condition describes 
a triangle in the complex plane whose area is 
proportional to the amount of CP violation in 
the model Following the observation of CP 
violation in the B° system [3], the focus has 
turned to testing the unitarity of the theory by 
over-constraining the sides and angles of this 
triangle. Most related measurements involve loop 
or box diagrams, and for which the CKM model is 
typically assumed when interpreting data [3] ■ This 
means the least- well determined observable, the 
phase 7 = arg (-V ud V* b /V cd V* b ) is of particular 
interest as 7 ^ can produce direct CP violation 
in tree decays. 

One of the most powerful methods for deter- 
mining 7 is measurements of the partial widths 
of B^ — > DK^ decays where the D signifies a 
D° or D° meson. In this case, the amplitude for 
the B~ —> D°K~ contribution is proportional to 
V c b whilst the B~ —> D°K~ amplitude depends 
on V u b- If the D final state is accessible for both 
D° and D° mesons, the interference of these two 
processes gives sensitivity to 7 and may exhibit 
direct CP violation. This feature of open-charm 
B~ decays was first recognised in its application to 
CP eigenstates, such as D — > K + K~ , 7r + 7r~ [5] but 
can be extended to other decays, e.g. D — > tt~K + . 
This second category, labelled "ADS" modes 
in reference to the authors of [BJ, requires the 
favoured, b — >• c decay to be followed by a doubly 
Cabibbo-suppressed D decay, and the suppressed 
b — >• u decay to precede a favoured D decay. The 
amplitudes of such combinations are of similar 
total magnitude and hence large interference can 
occur. For both the CP-mode and ADS methods, 
the interesting observables are partial widths and 
CP asymmetries. 

In this paper, we present measurements of 
the B^ decays in the CP modes, [K + K^}dIi ± 
and [iT + 'K~]]jh ± . the suppressed ADS mode 
[n^K^jjh^ and the favoured [K^tt^d^ combi- 
nation where h indicates either pion or kaon. De- 



cays where the bachelor — the charged hadron from 
the B~ decay — is a kaon carry greater sensitivity 
to 7. B~ — > Dtt~ decays have some limited sensi- 
tivity and provide a high-statistics control sample 
from which probability density functions (PDFs) 
are shaped. In total, 13 observables are measured: 
three ratios of partial widths 

R f T(B~ -» [f] D K-)+T(B+ -» [f] D K+) 
K l* T(B- [f] D 7r-) + T(B+ -+ [f] m +) ' 

(1) 

where / represents KK, tttt and the favoured Kit 
mode, six CP asymmetries 

,/ = r(s- -> [f]ph-) - T(B+ -» [f] D h+) 
h T(B- -> [f} D h~) + T(B+ -> [f] D h+) ' 1 1 

and four charge-separated partial widths of the 
ADS mode relative to the favoured mode 

h r(B± -> [K±nT] D h±Y { ) 

Elsewhere, similar analyses have established the 
B ± — > Dcph*- modes [7J 13 [9] and found evidence 
of the B ± -> [^K^d^ decay [lOl HH [Hj. 
Analyses of B ± -> [K^h+h^nK* decays [HE] 
have yielded the most precise measurements of 
7 though a 5cr observation of CP violation from 
a single analysis has not been achieved. This 
work represents the first simultaneous analysis of 
B ± -> Dcph^ and B ± — s- Z?ads'i ± modes. It is 
motivated by the future extraction of 7 which, 
with this combination, may be determined with 
minimal ambiguity. 

This paper describes an analysis of 1.0 fb _1 
of y/s = 7 TcV data collected by LHCb in 2011. 
The 2010 sample of 35 pb _1 is used to define 
the selection criteria in an unbiased manner. The 
LHCb experiment [15] takes advantage of the 
high bb and cc cross sections at the Large Hadron 
Collider to record large samples of heavy hadron 
decays. It instruments the pseudorapidity range 
2 < rj < 5 of the proton-proton (pp) collisions 
with a dipole magnet and a tracking system which 
achieves a momentum resolution of 0.4 — 0.6% in 
the range 5—100 GeV/c. The dipole magnet can 
be operated in either polarity and this feature is 
used to reduce systematic effects due to detector 
asymmetries. In 2011, 58% of data were taken 
with one polarity, 42% with the other. The pp 
collisions take place inside a silicon microstrip 
vertex detector that provides clear separation of 
secondary B vertices from the primary collision 
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vertex (PV) as well as discrimination for tertiary 
D vertices. Two ring-imaging Chcrcnkov (RICH) 
detectors with three radiators (aerogel, C4F10 
and CF4) provide dedicated particle identification 
(PID) which is critical for the separation of 
B~ — > DK~ and B~ — >• Dir~ decays. 

A two-stage trigger is employed. First a 
hardware-based decision is taken at a frequency 
up to 40 MHz. It accepts high transverse energy 
clusters in either an electromagnetic calorimeter or 
hadron calorimeter, or a muon of high transverse 
momentum (pt)- For this analysis, it is required 
that one of the three tracks forming the -B candi- 
date points at a deposit in the hadron calorimeter, 
or that the hardware-trigger decision was taken 
independently of these tracks. A second trigger 
level, implemented entirely in software, receives 

1 MHz of events and retains ~ 0.3% of them. It 
searches for a track with large px and large impact 
parameter (IP) with respect to the PV. This 
track is then required to be part of a secondary 
vertex with a high px sum, significantly displaced 
from the PV. In order to maximise efficiency at 
an acceptable trigger rate, the displaced vertex 
is selected with a decision tree algorithm that 
uses pt, Xipj flight distance and track separation 
information. Full event reconstruction occurs 
offline, and after preselection around 2.5 x 10 5 
events are available for final analysis. 

Approximately one million simulated events for 
each B ± — > [h + h~]Dh ± signal mode are used as 
well as a large inclusive sample of generic B — > 
DX decays. These samples are generated using a 
tuned version of Pythia [16] to model the pp col- 
lisions, EvtGen [T7] encodes the particle decays 
and Geant4 [TB] to describe interactions in the 
detector. Although the shapes of the signal peaks 
are determined directly on data, the inclusive sam- 
ple assists in the understanding of the background. 
The signal samples are used to estimate the rela- 
tive efficiency in the detection of modes that differ 
only by the bachelor track flavour. 

2 Event selection 

During event reconstruction, 16 combinations of 
B ± — > Dh^, D — > ft, ± ft, =F are formed with the can- 
didate D mass within 1765 - 1965 MeV/c 2 . 
D daughter tracks arc required to have 
Pt > 250 McV/c but this requirement is tightened 
to 0.5 < p T < 10 GeV/c and 5 < p < 100 GeV/c 



for bachelor tracks to ensure best pion versus kaon 
discrimination. The decay chain is refitted |19j 
constraining the vertices to points in space and 
the D candidate to its nominal mass, »ip DG [20]. 

Reconstructed candidates are selected using a 
boosted decision tree (BDT) discriminator |21) . 
It is trained using a simulated sample of B^ — > 
[K ± tt t ]dK ± and background events from the D 



sideband (35 < \m(hh) 



PDG 



< 100 MeV/c 2 



of the independent sample collected in 2010. The 
BDT uses the following properties of the candidate 
-B decay: 

• From the tracks, the D and B ± : px and xjp 
with respect to the PV; 

• From the B^ and D: decay time, flight dis- 
tance from the PV and vertex quality; 

• From the B^ : the angle between the momen- 
tum vector and a line connecting the PV to 
its decay vertex. 

Information from the rest of the event is employed 
via an isolation variable that considers the imbal- 
ance of Pt around the B^ candidate, 



ApT ~ MB) 



Pt(B) ~ E„PT 



(4) 



where the ^2 n PT sums over the n tracks within 
a cone around the candidate excluding the three 
signal tracks. The cone is defined by a circle 
of radius 1.5 in the plane of pscudorapidity and 
azimuthal angle (measured in radians). As no 
PID information is used as part of the BDT, it 
performs equally well for all modes considered here. 

The optimal cut value on the BDT response 
is chosen by considering the combinatorial back- 
ground level (b) in the invariant mass distribution 
of favoured S ± — > [Kt^dit^ candidates. The large 
signal peak in this sample is scaled to the antici- 
pated ADS-mode branching fraction to provide a 
signal estimate (s). The quantity s/V s + b serves 
as an optimisation metric. The BDT response 
peaks towards for background and 1 for signal. 
The optimal cut is found to be > 0.92 for the ADS 
mode; this is also applied to the favoured mode. 
For the cleaner CP modes, a cut of BDT > 0.80 
gives a similar background level but with a 20% 
higher signal efficiency. 

PID information is quantified as differences 
between the logarithm of likelihoods, ln£/,,, under 



2 



five mass hypotheses, h € {tt, K,p, e, fi} (DLL). 
Daughter kaons of the D meson arc required to 
have DLL^tt = In Ck — ln£ w > 2 and daughter 
pion must have DLL^ < — 2. Multiple candidates 
are arbitrated by choosing the candidate with the 
best-quality B+- vertex; only 26 events in the final 
sample of 157927 require this consideration. 

Candidates from B decays that do not contain 
a true D meson can be reduced by requiring the 
flight distance significance of the D candidate 
from the B~ vertex to be > 2. The effectiveness 
of this cut is monitored in the D sideband where 
it is seen to remove significant structures peaking 
near the B~ mass. A simulation study of the 
B~ -> K~K + K~, K~ir + ir~ and K~ K+ tt~ modes 
suggests this cut leaves 2.5, 1.3 and 0.8 events 
respectively under the B~ — > [K + K ~]dK ~, 
[-k + -k~]dK~ and [it + K~]dK~ signals. This cut 
also removes cross feed {e.g. B~ — > [K ^tt + ]d^ 
as a background of [ti + 'k^]£)K~) which occurs 
when the bachelor is confused with a D daughter 
at low decay time. Finally, the combination of 
the bachelor and the opposite-sign D° daughter is 
made under the hypothesis they are muons. The 
parent B candidate is vetoed if the invariant mass 
of this combination is within ±22 MeV/c 2 of either 
the J/ip or i/)(2S) mass [21]] . 

Due to misalignment, the reconstructed B+- 
mass is not identical to the established value, 



m 



PDG 



|. As simulation is used to define 
background shapes, it is useful to apply linear 
momentum scaling factors separately to the two 
polarity datascts so the B+- mass peak is closer 
to Tip DG . After this correction, the D° — > K~ir + 
mass peak is measured at 1864.8 MeV/c 2 with a 
resolution of 7.4 MeV/c 2 . Selected D candidates 
are required to be within ±25 MeV/c 2 of mp DG . 
This cut is tight enough that no cross feed occurs 
from the favoured mode into the CP modes. In 
contrast, the ADS mode suffers a potentially large 
cross feed from the favoured mode in the circum- 
stance that both D daughters arc misidentified. 
The invariant mass spectrum of such cross feed 
is broad but peaks around nip DG . It is reduced 
by vetoing any ADS candidate whose D candidate 
mass under the exchange of its daughter track mass 
hypotheses, lies within ±15 MeV/c 2 of mp DG . 
Importantly for the measurements of R^, this 
veto is also applied to the favoured mode. With 
the D mass selection and the D daughter PID 
requirements, this veto reduces the rate of cross 



feed to an almost negligible rate of (6 ± 3) x 10 5 . 

Partially reconstructed events populate the in- 
variant mass region below the B+- mass. Such 
events may enter the signal region, especially where 
Cabibbo-favoured B — > XDtt^ modes are misiden- 
tified as B — > XDK+-. The large simulated sam- 
ple of inclusive B q — > DX decays, q € {u, d, s}, 
is used to model this background. After applying 
the selection, two non-parametric PDFs [35] are de- 
fined (for the Dtt^ and DK+ selections) and used 
in the signal extraction fit. These PDFs are ap- 
plied to all four D modes though two additional 
contributions are needed in specific cases. In the 
D -t K+K- mode, A° b -> [p+ A'"7r+] Ac ^ enters if 
the pion is missed and the proton is reconstructed 
as a kaon. In the B ± — > DadsK^ mode, partially 
reconstructed B® — > D K + tt~ decays represent an 
important, Cabibbo-favoured background. PDFs 
of both these sources are defined from simulation, 
smeared by the modest degradation in resolution 
observed in data. When discussing these contribu- 
tions, inclusion of the charge conjugate process is 
implied throughout. 

3 Signal yield determination 

The observablcs of interest are determined with 
a binned maximum-likelihood fit to the invariant 
mass distributions of selected B candidates [23] . 
Sensitivity to CP asymmetries is achieved by sep- 
arating the candidates into B~ and B + samples. 
B ± — > DK+- events are distinguished from B+ — > 
using a PID cut on the DLLk-k of the bachelor 
track. Events passing this cut are reconstructed as 
DK+, events failing the cut are reconstructed as 
the final state. The fit therefore comprises 

four subsamples — (plus, minus) x (pass, fail) — for 
each D mode, fitted simultaneously and displayed 
in Figs. HHU The total PDF is built from four or 
five components representing the various sources of 
events in each subsample. 

1. B+ — > . In the sample failing the bach- 
elor PID cut, a modified Gaussian function, 



f(x) oc exp 



2cr 2 + (x- n) 2 a L , R 



(5) 



describes the asymmetric peak of mean \x and 
width a where ai(a; < fi) and olr(x > fi) pa- 
rameterise the tails. 

True B± -t _D?r ± events that pass the PID 
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cut arc reconstructed as B ± — s- DK^. As 
these events have an incorrect mass assign- 
ment they form a displaced mass peak with 
a tail that extends to higher invariant mass. 
These events are modelled by the sum of two 
Gaussian PDFs also altered to include tail 
components. All parameters are allowed to 
vary except the lower-mass tail which is fixed 
to ensure fit stability and later considered 
amongst the systematic uncertainties. These 
shapes are considered identical for B~ and 
B + decays and for all four D modes. This 
assumption is validated with simulation. 

2. B^ — > DK^ : In the sample that passes the 
DLL^ cut on the bachelor, the same modi- 
fied Gaussian function is used. The mean and 
the two tail parameters are identical to those 
of the larger, B ± — >• Dtt^ 1 peak. The width 
is 0.95 ±0.02 times the Dn^ width, as deter- 
mined by a standalone study of the favoured 
mode. Its applicability to the CP modes is 
checked with simulation and a 1% system- 
atic uncertainty assigned. Events failing the 
PID cut are described by a fixed shape that 
is obtained from simulation and later varied 
to assess the systematic error. 

3. Partially reconstructed B —> DX: A fixed, 
non-parametric PDF, derived from simula- 
tion, is used for all subsamples. The yield in 
each subsample varies independently, making 
no assumption of CP symmetry. 

4. Combinatoric background: A linear approx- 
imation is adequate to describe the slope 
across the invariant mass spectrum consid- 
ered. A common parameter is used in all sub- 
samples, though yields vary independently. 

5. Mode-specific backgrounds: In the D — > KK 
mode, two extra components are used to 
model A° — > A+/i~ decays. Though the to- 
tal contribution is allowed to vary, the shape 
and relative proportion of A~£K~ and A+7r~ 
arc fixed. This latter quantity is estimated at 
0.060 ±0.01 5, similar to the effective Cabibbo 
suppression observed in B mesons. For the 
B ± — > Dads-^ ± mode, the shape of the 
B® — > D°K + ir~ background is taken from 
simulation. In the fit, this yield is allowed 
to vary though the reported yield is consis- 
tent with the simulated expectation, as de- 
rived from the branching fraction |24) and the 
bb hadronisation [25] . 



The proportion of B — > Dh^ passing or 
failing the PID requirement is determined from a 
calibration analysis of a large sample of D*^ de- 
cays reconstructed as D**- — >• Dtt^, D — s- K^^. 
In this calibration sample, the K and 7r tracks may 
be identified, with high purity, using only kine- 
matic variables. This facilitates a measurement 
of the RICH-based PID efficiency as a function 
of track momentum, pseudorapidity and number 
of tracks in the detector. By reweighting the 
calibration spectra in these variables to match 
the events in the _B ± — > Dn^ peak, the effective 
PID efficiency of the signal is deduced. This 
data-driven technique finds a retention rate, for a 
cut of DLLktt > 4 on the bachelor track, of 87.6% 
and 3.8% for kaons and pions, respectively. A 
1.0% systematic uncertainty on the kaon efficiency 
is estimated from simulation. The B ± — > Z?7r ± fit 
to data becomes visibly incorrect with variations 
to the fixed PID efficiency > ±0.2% so this value 
is taken as the systematic uncertainly for pions. 

A small negative asymmetry is expected in 
the detection of A~ and A' + mesons due to their 
different interaction lengths. A fixed value of 
(— 0.5 ± 0.7)% is assigned for each occurrence of 
strangeness in the final state. The equivalent 
asymmetry for pions is expected to be much 
smaller and (0.0 ± 0.7)% is assigned. This un- 
certainty also accounts for the residual physical 
asymmetry between the left and right sides of 
the detector after summing both magnet-polarity 
datasets. Simulation of B meson production in pp 
collisions suggests a small excess of B + over B~ 
mesons. A production asymmetry of (—0.8 ±0.7)% 
is assumed in the fit such that the combination 
of these estimates aligns with the observed raw 
asymmetry of B ± — > J/ipK^ decays at LHCb |26j . 
Ongoing studies of these instrumentation asym- 
metries will reduce the associated systematic 
uncertainty in future analyses. 
The final B ± — > Dh^ 1 signal yields, after summing 
the events that pass and fail the bachelor PID cut, 
are shown in Table [TJ The invariant mass spectra 
of all 16 B ± -> [h+h-] D h ± modes are shown 
in Figs. HHU Regarding the B ± —> Z?7r ± mass 
resolution: respectively, 14.1 ± 0.1, 14.2 ± 0.1 and 
14.2 ± 0.2 MeV/c 2 are found for the D -> KK, 
Kn and tttt modes with common tail parameters 
a L = 0.115 ± 0.003 and a R = 0.083 ± 0.002. As 
explained above, the B ± — > DK^ widths are fixed 
relative to these values. 
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Table 1: Corrected event yields. 



B mode D mode 


B~ B+ 


DK ± 


3170 ± 83 3142 ± 83 
592 ± 40 439 ± 30 
180 ± 22 137 ± 16 
23 ± 7 73 ± 11 


JJtt ± T 

TV 7T T 
TV+K* 


40767 ± 310 40774 ± 310 
6539 ± 129 6804 ± 135 
1969 ± 69 1973 ± 69 
191 ± 16 143 ± 14 



The ratio of partial widths relates to the ratio 
of event yields by the relative efficiency with which 
£? ± — > D°K ± and B^ — > D°7r ± decays are recon- 
structed. This ratio, estimated from simulation, 
is 1.012, 1.009 and 1.005 for D KK, Kit, hit 
respectively. A 1.1% systematic uncertainty, based 
on the finite size of the simulated sample, accounts 
for the imperfect modelling of the relative pion 
and kaon absorption in the tracking material. 

The fit is constructed such that the observables 
of interest are parameters of the fit and all sys- 
tematic uncertainties discussed above enter the fit 
as constant numbers in the model. To evaluate 
the effect of these systematic uncertainties, the fit 
is rerun many times varying each of the system- 
atic constants by its uncertainty. The resulting 
spread (RMS) in the value of each observable is 
taken as the systematic uncertainty on that quan- 
tity and is summarised in Table [2j Correlations 
between the uncertainties are considered negligible 
so the total systematic uncertainty is just the sum 
in quadrature. For the ratios of partial widths in 
the favoured and CP modes, the uncertainties on 
the PID efficiency and the relative width of the 
DK ± and Dtt^ 1 peaks dominate. These sources 
also contribute in the ADS modes, though the as- 
sumed shape of the B° s — > D°K + ir~ background 
is the largest source of systematic uncertainty in 
the B ± j D A ds^ ± case. For the CP asymme- 
tries, instrumentation asymmetries at LHCb are 
the largest source of uncertainty. 



Table 2: Systematic uncertainties on the observ- 
ables. PID refers to the fixed efficiency of the 
DLLa'tt cut on the bachelor track. PDFs refers to 
the variations of the fixed shapes in the fit. "Sim" 
refers to the use of simulation to estimate relative 
efficiencies of the signal modes which includes the 
branching fraction estimates of the A° background, 
^instr. quantifies the uncertainty on the production, 
interaction and detection asymmetries. 



xl0~ 3 


PID 


PDFs 


Sim 


-^instr. 


Total 


r>Kw 

U K/tt 


1.4 


0.9 


0.8 





1.8 


r/KK 

U K/tt 


1.3 


0.8 


0.9 





1.8 


T? nn 
U K/w 
A Kit 


1.3 


0.6 


0.8 





1.7 





1.0 





9.4 


9.5 




0.2 


4.1 





16.9 


17.4 


aKK 
A K 


1.6 


1.3 


0.5 


9.5 


9.7 


A nn 
A K 


1.9 


2.3 





9.0 


9.5 


aKK 


0.1 


6.6 





9.5 


11.6 


ATTTT 

^ir 


0.1 


0.4 





9.9 


9.9 


R K 


0.2 


0,1 





0.1 


0.4 


R K 


0.4 


0.5 





0.1 


0.7 


R TT 


0.01 


0.03 





0.07 


0.08 


Rt 


0.01 


0.03 





0.07 


0.07 



4 Results 

The results of the fit with their statistical uncer- 
tainties and assigned systematic uncertainties are: 



n K/TT — 





0774 ±0.0012 ±0.0018 


U K/tv — 





0773 ±0.0030 ±0.0018 


U K/ir — 





0803 ±0.0056 ±0.0017 


A?" = 


-0 


0001 ±0.0036 ±0.0095 


A«* = 





0044 ±0.0144 ±0.0174 


aKK 

A K — 





148 ± 0.037 ±0.010 


A IC = 





135 ±0.066 ±0.010 


A« K = 


-0 


020 ±0.009 ±0.012 


A7 = 


-0 


001 ±0.017 ±0.010 


R~ K = 





0073 ± 0.0023 ±0.0004 


R+ K = 





0232 ±0.0034 ±0.0007 


K = 





00469 ±0.00038 ±0.00008 


Rt = 





00352 ± 0.00033 ± 0.00007 



From these measurements, the following quantities 



■5 



can be deduced: 
Rcp+ ~ 



.4 



CP+ 



< n K/wi U K/Tr > / U K/tv 

1.007 ± 0.038 ±0.012 

7T7T 



<A« K ,A 



K 



> 



R 



•ADS(K) 



A 



ADS(K) 



-RaDS(tt) 



.4 



ADS(tt) 



0.145 ± 0.032 ±0.010 
(R K +R+)/2 

0.0152 ± 0.0020 ±0.0004 
(R K -R+)/(R K + R+) 
-0.52 ±0.15 ±0.02 
(R- + R+)/2 

0.00410 ± 0.00025 ± 0.00005 
(R- - R+)/(R~ + R+) 

0.143 ± 0.062 ±0.011, 



where the correlations between systematic uncer- 
tainties are taken into account in the combination 
and angled brackets indicate weighted averages. 
The above definition of Rcp+ is only approximate 
and is used for experimental convenience. It as- 
sumes the absence of CP violation in B^ — > Dtt^ 
and the favoured B ± — > DK ± modes. The exact 
definition of Rcp+ is 



T(B- Dcp+K-) + T(B+ Dcp+K^ 



T(B- 



D°K~) 



(0) 



so an additional, and dominant, 1% systematic un- 
certainty accounts for the approximation. For the 
same reason, a small addition to the systematic un- 
certainty of R^Jn is needed to quote this result as 
the ratio of B^ branching fractions, 



B(B~ 



D°K~) 

- — f = (7.74 ± 0.12 ± 0.19)%. 



B(B~ D°tt- _ 

To summarise, the B ± — > DK^ ADS mode is 
observed with ~ 10cr statistical significance when 
comparing the maximum likelihood to that of 
the null hypothesis. This mode displays evidence 
(4. Oct) of a large negative asymmetry, consistent 
with the asymmetries reported by previous exper- 
iments PUCEUCEi]. The B ± -» L»tt ± ADS mode 
shows a hint of a positive asymmetry with 2.4er 
significance. The KK and ttit modes both show 
positive asymmetries. The statistical significance 
of the combined asymmetry, Acp+, is 4.5<r which 
is similar to that reported in |U] albeit with a 
smaller central value. All these results contain 
dependence on the weak phase 7 and will form an 
important contribution to a future measurement 
of this parameter. 



Assuming the CP-violating effects in the CP 
and ADS modes are due to the same phenomenon 
(namely the interference of b — > cus and b — > ucs 
transitions) we compare the maximum likelihood 
with that under the null-hypothesis in all three D 
final states where the bachelor is a kaon. This log- 
likelihood difference is diluted by the non-negligible 
systematic uncertainties in Acp+ and Aads(.R") 
which are dominated by the instrumentation asym- 
metries and hence are highly correlated. In conclu- 
sion, with a total significance of 5.8<r, direct CP 
violation in B ± — s- DK ± decays is observed. 

Acknowledgements 

We express our gratitude to our colleagues in 
the CERN accelerator departments for the ex- 
cellent performance of the LHC. We thank the 
technical and administrative staff at CERN and 
at the LHCb institutes, and acknowledge sup- 
port from the National Agencies: CAPES, CNPq, 
FAPERJ and FINEP (Brazil); CERN; NSFC 
(China); CNRS/IN2P3 (France); BMBF, DFG, 
HGF and MPG (Germany); SFI (Ireland); INFN 
(Italy); FOM and NWO (The Netherlands); SCSR 
(Poland); ANCS (Romania); MinES of Russia and 
Rosatom (Russia); MICINN, XuntaGal and GEN- 
CAT (Spain); SNSF and SER (Switzerland); NAS 
Ukraine (Ukraine); STFC (United Kingdom); NSF 
(USA). We also acknowledge the support received 
from the ERC under FP7 and the Region Au- 
vergne. 



G 




m{D}f) (MeV/c 2 ) 

Figure 1: Invariant mass distributions of selected S ± — > [if ± 7r =F ]i)/i ± candidates. The left plots are B~ 
candidates, B + are on the right. In the top plots, the bachelor track passes the DLL/^ > 4 cut and the B 
candidates are reconstructed assigning this track the kaon mass. The remaining events are placed in the 
sample displayed on the bottom row and are reconstructed with a pion mass hypothesis. The dark (red) 
curve represents the B — > DK^ 1 events, the light (green) curve is B — > D^. The shaded contribution 
are partially reconstructed events and the total PDF includes the combinatorial component. 
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Figure 2: Invariant mass distributions of selected B ± — > [K + K }oh ± candidates. See the caption of 
Fig. [T]for a full description. The contribution from A& — > Afh T decays is indicated by the dashed line. 
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m{Dhr) (MeV/c 2 ) 

Figure 4: Invariant mass distributions of selected -B ± — >■ [7r ± iC T ]i3/i ± candidates. See the caption of 
Fig. Q] for a full description. The dashed line here represents the partially reconstructed, but Cabibbo 
favoured, B® — > D°K~ir + and — ^ D°K + ir~ decays where the pions are lost. The pollution from 
favoured mode cross feed is drawn, but is too small to be seen. 
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